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I. Introduction 


This paper presents an analytical model for a free convective 
boundary layer diffusion flame and compares results from this mode! to 
experiment. The analytical model will add two new features to the solution 
presented by Pagnil1] in order to improve agreement with experimental 
data [2]. The analysis builds upon the work originally done by Kosdon, 
Williams, and Buman [3]; and Kim, de Ris, and Krosser [4]. It is proposed 
that the fuel undergoes endothermic decomposition before reaching the 
flame. Several experimental studies have demonstrated this phenomena 
[ 5,6,7]. The endothermic decomposition of the fuei will be referred to as 
cracking. The effect of cracking on the free convective diffusion flame is 
explored by the use of modified Shvab Zeldovich variables which mode! both 
the flame reaction and the cracking reaction [8]. The second addition to the 
solution given by Pagni is the use of a temperature dependent density - 
viscosity (pn) function [9] which allows the ambient viscosity to be 
determined by ambient conditions. In order to match the location of the 
theoretical flame with the location of the experimental flame, previous 
studies [10] have adjusted the ambient viscosity to non ambient values. By 
using the temperature dependent density-viscosity (pp) function, the 
analytical model accurately predicts the location of the flame using the 


ambient value of viscosity. 


Figure | shows a schematic diagram of the free convective diffusion 
flame. The wall is a pyrolyzing slab from which fuel is endother mically 


pyrolyzed by heat flux to the wall. The fuel then diffuses toward the flame. 


At a sufficiently high temperature the fuel undergoes an endothermic 
reaction in which the fuel molecules Aecontore into simpler molecules. This 
reaction zone will be modeled as a sheet in order to simplify the analysis. 
The actual cracking reaction occurs over a fairly wide region. The cracked 
fuel then diffuses toward the flame and toward the wall. Cracked fuel 
‘reaching the flame burns in an exothermic reaction with orygen generating 
products and providing heat to drive both the flow and the combustion 


process. 
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Fig. 1. 


The analysis of the flame incorporating the effects of thermal cracking 
and temperature dependent density viscosity product (pp) occupies the first 
half of the paper. The analysis begins with the governing equations and 
their boundary conditions. Following this, the assumptions used to simplifv 


the governing equations are given. The analysis proceeds with the 


B, 
derivation of the modified Shvab Zeldovich variables used to eliminate the 
generation terms in the energy and species equations. The determination of 
enthalpy and species fields necessary for the solution of the problem are 
presented next. Nondimensional groups describing the flow field are then 
introduced which include a Howarth-Dorodnitzyn transformation to make 
the governing equations appear incompressible. The variable properties are 
introduced at this point. The problem is reduced from two dimensions to 
one dimension through the use of a similar stream function. The governing 
equations are then restated in their simplified form. Siz nondimensional 


groups are shown to describe the system. 


The second half of the paper examines the results of the analytical 
model. Solutions with and without cracking are compared with 
experimental data [2] for temperature and streamwise velocity. The fuels 
modeled are cyclohexane, ethanol, iso-octane, and n-heptane. The pyrolysis 
rate at the wall is predicted for each of the fuels based on the theoretical 
results. Finally, the nondimensional groups describing the cracking reaction 
are varied in a parametric study to determine there effect on the flame 


solution. 


I]. Analysis 
GC ee Raval 


Symbols are defined in the nomenclature. The flame is assumed to be 
steady and two dimensional with coordinate system shown in figure 1. With 


the addition of fuel cracking, the one step reaction, 


Pear sumer eres {36 MESS sin BiPracs aEkhy Hae & areD, 


F/ 


is broken up into two reactions. The first reaction is the cracking reaction : 
i) 
Ve, ie , seks Ve > re - ier . (2) 
And the second reaction is the modified flame reaction in which the cracked 


fuel reacts with the oridizer: 


ybiites + Com (ol eicten Ohanete alle 


The enthalpies of reaction are related by 


Qerec 7 Cs. us Q er é (4) 


In the analysis only Qc- and Qr will be used. Where Qcrr; is needed, the 


right hand side of eqn. (4) will be used. 


The conservation equations and boundary conditions are as follows. 


Conservation of mass: 


Vp aoe (5) 
Conservation of momentum : 

oa. Vee ~ Winter) = gif) = ye 
with boundary conditions for streamwise velocity given by 

Una oa U Cony oat= ANC Sy oe eee ee b7) 


The boundary conditions for transverse velocity account for the blowing at 


the wall. vide 3 
WV (%5 9) = — and V7 (0,4) = ANC, ed Fey 


(8) 


Conservation of energy : 


ee fey 2 acres) deeh al Tee vragen linw 907) 


where 


a. 


N) 


‘ 
lo Sy gow | (/0) 


The energy equation contains two source terms. The first term, GFL. 


3! 


represents the heat produced by the flame reaction. The second source term. 


ace accounts for the energy absorbed by the cracking reaction. The 


boundary conditions for the energy equation are : | 
| he neve and h (0,;y) = EE OTe a 0. (WW 


The species equation must be solved for the uncracked fuel ( Yr}), the 
cracked fuel ( Yr), and the oxidizer ( Y,). The species equation for Yq is 


path — (20K) = my (12) 


with boundary conditions 


pu(Gere) =o and YiCo,y) = Le OS leone 


The uncracked fuel species equation is 


as 


saan he be —U+(20%) een ae 


where mrjcp represents the mass sink of uncracked fuel at the cracking 


sheet. Ambient boundary conditions for the uncracked fuel are 
homogeneous. The wall! boundary condition is determined by a mass balance 


at the wall where Yrry is used to account for the inerts in the pyrolyzed 


material. j O Va, 


rn. ee 5] mM, Le i> 20 In ly ‘ lis) 


The cracked fuel species equation is 


par. Vie V- (2D Me ) a Cosy 1 SOUT (16) 
where myocp is a source of cracked fuel at the cracking reaction, and mppy is 
a sink of cracked fuel at the flame reaction. Ambient boundary conditions 
for the cracked fuel are also homogeneous. The wall boundary condition for 


cracked fuel assumes that the cracked fuel does not diffuse into the wall. 


d Ver 
Pe Elune eae 


(17) 
An energy balance at the wall determines the wall mass flux. 
js eee 
w L *h a an o cae 


This completes the list of governing equations. The introduction of the 
variable properties will be introduced when the equations are 


nondimensionalized. 


NI 


3B. Assumptions 


The following are a list of simplifications and assumptions used tn the 


analytic model. 


1. The flow is assumed to be steady, laminar, and two dimensional. 


2. Diffusion is governed by Fick's Law with the same diffusion 
coefficients for all species. 


2. The overall reaction will be modeled as a two step reaction: a 
flame reaction and a cracking reaction occurring at two separate 
regions in the boundary layer. 


3. The flame reaction zone and the cracking reaction zone will be 
modeled as sheets. 


4. All undecomposed fuel which reaches the cracking sheet will be . 
assumed to decompose fully and irreversibly into cracked fuel. 


5. The cracked fuel cannot diffuse into the pyrolyzing slab. 


6. Cracked fuel and oxygen are completely consumed in stoichiometric 
proportions at the flame sheet to form products. 


7. The cracking reaction takes place at a known temperature. 
8. The Lewis number equals one. 
9. Specific heat is not a function of temperature. 


10. The density viscosity product is modeled using the viscosity of 
Nitrogen. 


11. The wall temperature is the boiling temperature of the fuel at 
atmospheric pressure. 


12. Assume a similar solution. 
13. Neglect radiation. 


14. Ambient is at rest. 


The first step in the analysis is the removal of the generation terms in 
the energy and species equations. By assuming unity Lewis number and 
stoichiometric consumption of species, it 1s possible to eliminate the 
generation terms by using modified Shvab Zeldovich variables. The 
treatment of multiple flame sheets by the use of Shvab Zeldovich variables 
was originally developed by Penner and Libby [8]. Equation (2) limits the 


consumption of energy and species at the cracking sheet to 


at oo” 


- My = Pr arcr — para hh ¢ 7) 
Ve Me Ve AP Qer , 


Similarly, Equation (3) limits the consumption of energy and species at the 


eo 7 PRL/A = 
flame sheet to _17,,, ii Deca ee oun fire (29) 


Ve “Nr i? Yo Mo ee Ger 


By using eqns. (19) and (20), the species equation for Yo, Yr, and Yro can be 


reduced to a single, homogeneous equation with the independent variable 
Yorn th pw welt» gedos ueike 


los} = VeUe Vp Me Vounto 


(2) 


This variable is normalized such that the boundary conditions are either 


unity or zero by b- Gv 


= on 
Tar—a og eae E ris 


The species mass fractions are nondimensionalized by the ambient mass 
fraction of oxygen and the stoichiometric ratio 
Perot We Pr ge Niet tae V2 ve MN and ys (23) 
hs Youu Ve lr “a Yoro Ve Mp Yo 2 


The nondimensional uncracked fuel at the wall will be referred to as the 


uncracked fuel equivalence ratio 
Yai: Yo A, 


p = genpere ; | (29 
Fi 1 Ve rhe 


The cracked fuel at the wall will be refered to as the cracked fuel 
equivalence ratio 

be 

g a Mike BOGS Case 

Fa A pe Ve (Nn - 

These variables are unknown at this point. Using the nondimensiona! 
variables, the normalized Shvab Zeldovich variable becomes 

qisv2 SATA. 3 | (26) 

Pei % Pe,’ ‘if / , 

and the species equation reduces to 


pre a (oP Vp O Q7) 


with boundary conditions 
Tryp Od raF ord ane aL Oy) = J (x,2) = Me is (28) 


The generation term in the energy equation is eliminated using the following 


independent variable 


L pling veers a Mee (= es Pir oe 
yor Bae ag Ot MVM Tid Mh me vempne As 
This variable is normalized such that 
fa ARS (3o) 


Ss A eT ne ye 
With the specific heat assumed constant, the enthalpy 1s 


ba ee we orice mmemae oa (3) 


The ambient oxygen and the heat released by the overall flame reaction are 
used to nondimensionalize the enthalpy so that the dimensionless enthalpy 
is defined as = ot Sepia pa Yow Line, 

SA eefefe (Se 
Q ee Yo 2 < 


ae 
| 


This variable is of order one at the flame. The dimensionless enthalpy at the 


wall is defined using eqn.(32)as ~ _ Cp (72 - Ta) vw ™, 
hw a @ e A 
Fe flay 


The ratio of the cracking enthalpy of reaction to the overall enthalpy 
of reaction indicates the relative magnitude of the heat absorbed by the 
cracking reaction and is defined as 


~ Goer 
= — 5 
Cre i ( y) 


Using the previously defined nondimensional groups, the normalized energy 


species variable as defined by eqn. (30) becomes 


Tt EDS (35) 
hy oo EN 


The differential equation for the energy species variable is eqn. (27) with 


boundary conditions given by eqn. (28). 


In the solution to the problem with cracking, the parameters used to 
normalize the species and the energy species variables ( Yrj~ .Y¢oy,) are 
variables. In the solution of the problem without cracking, these parameters 


could be determined 4 prvor/ by the mass transfer number,B , and the mass 


1] 


fraction of fuel in the transferred gas, Yry. The existence of the endothermic 


cracking sheet causes the mass fraction of uncracked fuel at the wall to be 
dependent upon both the mass fraction of cracked fuel at the wall and the 
heat absorbed by the cracking sheet. The traditional mass transfer number 
is defined as 


Boizon| fsa A oot hnaloa hinoloviur ony Cd 


and the nondimensional fuel in transferred medium is defined as 


dh peo = Yer Yo M. 


y ey eae: 
ogo - F 


Eor the problem without cracking, o¢; is determined by the energy and mass 
balance at the wall to be : 


ay 
or, an Wesus el ; (38) 


Bere ee ee 


With the addition of cracking, the mass transfer number is modified by a 
cracking mass transfer number defined as 
B 2 Yeow Ger (39) 
at Vamp (pe aes 


An effective mass transfer number is then defined as 


Bhined Rlieary = id 


The value of or; is determined by eqn (15), eqn (17), and eqn(18) for the 
energy and mass balances at the wall and eqn. (26) and eqn. (35) for the 


definitions of the Shvab Zeldovich variables giving 


d ri ait ee ae Cnr ete Pete ae > (+/) 


ep ae a, 


12 


The value of of is still unknown and must be determined by the solution of 


the cracked fuel species equation inside the cracking sheet. 


Desnec | Enthalpv Field 
In order to solve for the species and enthalpy fields, it is necessary to 
solve the differential equation for J. A system of iteration is needed since 
the solution for J is dependent upon the enthalpy field through the 
momentum equation. In this section, it is assumed that the solution for J is 
known. The following results rely on the assumption that all species 
diffusing toward either the flame sheet or cracking sheet are completely 
consumed. From the cracking sheet to the ambient environment there are 
only two field variables present: either the cracked fuel and enthalpy, or the 
oxidizer and enthalpy. Since there are two independent Shvab Zeldovich 
variables and only two field variables, the enthalpy and species fields are 
completely specified by the solution to J. Between the cracking sheet and 
the wall, both the cracked fuel and the uncracked fuel are present with the 
enthalpy. Since there are three field variables and only two Shvab Zeldovich 
Variables, it is necessary to solve the species equation for the cracked fuel in 


the region between the cracking sheet and the wall. 


The location of the cracking sheet is specified by the temperature at 
which the cracking reaction occurs. It is also know that the mass fraction of 
the uncracked fuel goes to zero at the cracking sheet. The location of the 
cracking sheet in the solution for J is found from this information. With the 


assumption of constant specific heat, the dimensionless enthalpy at the 


cracking sheet defined as 
ae Cp ( ioe ne ix ) A MY, 4 


h S ! Le 12.) 
ue rye Ves 


The location of the cracking sheet in J is determined from the species and 


the energy species variables as 
lo+ oe - ie 
SS ° (Pape) 
ck pee OO — en Q - he 


By solving the species variable for yr and using the location of the cracking 


sheet in J, it is then possible to determine the value of yr2,, at the cracking 


a 7 = Jer [ Peas Pr i ai | ° (+4) 


YW recrR 


The boundary conditions for the cracked fuel species equation are now 
completely specified by the zero flux condition at the wall, eqn.(17), and the 


value of Yroc, at the cracking sheet, eqn.(44). 


The location of the flame sheet in ] is found by setting y, and yp 


equal to zero so that eqn. (26) reduces to 
| (45) 
on te Pees 
The location of the flame sheet in J separates the region in which the oxidizer 
is zero from the region in which the cracked fuel is zero. With this — 
information it now possible to determine the enthalpy field. The enthalpy 
between the cracking sheet and the wall is determined by 


o> zh, Reeling anpRtrepdgneec i ea 


The enthalpy between the cracking sheet and the flame sheet is given by 
~ Si Oana bea neda ain Aa eee 


and enthalpy from the flame sheet to ambient is 
LY ey ane mnd ToL ceria Heche: (¥ 7) 
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E. Treatment of the Boundary Layer Equations 

The treatment of the boundary layer equations follows that given in 
previous work [2] except for the addition of a temperature dependent 
density viscosity product (pp). The density viscosity product is determined 
using the viscosity of nitrogen and is calculated using the.collision integrals 


for a Leonard-Jones potential given in Hirschfelder, Curtis and Bird [9]. 


Curve fitting density viscosity product with T / Ting yields 


r at AS PR 


Pr “co 


The boundary layer equations for ] and the cracked fuel are the same. 


For the cracked fuel equation, J is replaced by 3 


WR RUOBW Ie Tani 2o( eo Bale ee) 


The momentum equation is 
o ; 
ay = ‘pk jee al at a — 904 --@), (YD 


A Howarth-Dorodnitzyn transform 1s used to make the governing equations 
appear incompressible. The nondimensional variables used to describe the 


flow field are given as, 


ly va | 
x Lies ce ; CGa) 
5 7 A P: Ss % ‘ae cle ~G 7 
va Be 
=a and V = 5 —— (2+ 4 dy) , 
yi Gr,* Vio Or, fro “A fz 


V eyed 
- + a =O 5 Say 
y EE GRRE perdi oe: dhivate BEY 
ey ev ee Me ae ig nr aera le) 

A stream function is used to satisfy conservation of mass. 
AO 22 and PE Ge a ' (ro 


The solution is assumed to be similar by the use of a similar stream function 


lame aieemamhe); (5 7) 


and a similarity variable defined as 


Sex 
ALN aE EUR (52) 


This reduces the problem from two dimensions to only one dimension. - 


educed Governi ions 


The governing equations can now be restated in there reduced form. 


The differential equation for the energy species variable is given as 


Pic DLS pte: Po alae 8 told 


with boundary conditions 


Tees | ad» WW) Come ne (60) 


The cracked fuel species equation inside the cracking sheet is specified by 
rae et i b 
PCN | + Pou yore a 


with the zero flux into the wall condition given as 


Oo Ver 


8 Pr P(e) eee = Nee 


¢ 


A= 0 (62) 


and the boundary condition at the cracking sheet given by 


Yrr itasy ) a Y F2CR . (63) 


where eqn. (44) is used to determine y,, «, he value of ner is found by 
searching the solution of J(n ) for J,,as specified by eqn. (43). A similar 
method is used to find nrj. The reduced momentum equation is 


paris Master me Ey eet ca = ea ‘ete 


with boundary conditions 


/a2s Sein yton Da vag ees op IEC Ee (65) 
The blowing condition at the wall is given by 
3h oe) = bew oe (66) 


where Berr is defined by eqn(40). The enthalpy field is given by eqn. (46), 
(47) , and (48). | | 


G. § li CS ’ , 


The system of equations are highly coupled and must be solved in an 
iterative manner. The equations that must be solved are as follows: 
Equations (59) through (66) give the governing equations and main 
boundary conditions; eqns. (43) through (49) determine the enthalpy field, 
the density viscosity product and one of the boundary conditions for the 
cracked fuel equation: and eqns. (39) through (41) determine the effective 


mass transfer number and the uncracked fuel equivalence ratio. 


The diffusion flame with cracking is controlled by six nondimensional 
parameters listed in table 1. Besides these parameters, the Prandtl number 
and the ambient properties of air must be known. The effect of cracking 
adds two new parameters to the analysis: the heat of reaction ratio gives the 


ratio of the cracking heat of reaction to the total heat of reaction and 


[8 


indicates the relative amount of heat absorbed by the cracking sheet. The 
cracking enthalpy determines the location of the cracking sheet in the 


enthalpy field. 


The problem 1s solved numerically using the C programing language. 
The momentum equation is linearized and solved using a tridiagonal solver. 
The energy-species equation and cracked fuel equation are solved using 
senaration of variables and integration. The complex interdependencies are 
solved by iteration. 
TABLE 1 


Summa-ryof Dimensionless Parameters 


Fuel in Transferred medium 
= Lor vo Mle 
fy = Yeas tt 


Mass transfer number 


Qre Te xe ay 
ba Me. 4, Ae | a A. 
Wall enthalpy 
~— Cp ( ifs a ay Yo JY 5 
ip. eis anaes 


Cracking enthalpy 
G if > ae es ) vo he 


h cA Fi a ate lace 
Heat of reaction ratio 
~ = Dp 
i Gee 


Grashof number evaluated at 7 
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Paes Pe NG = 
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Discussion and Results 
A. Fuels Modeled 


The fuels chosen for analysis were based on the availability of 
experimental data for temperature. and velocity [2] fields. The fuels 
modeled were n-heptane, iso-octane, cyclohexane, and ethanol. The 


general overall reaction for these fuels can be written as 


Vik om a as NEN + plete ar ee (67) 


where the stoichiometric coefficients, the fuel molecules, and the 
enthalpies of reaction are given in table 2. At the cracking reaction, the 
fuels were assumed to decompose into acetylene ( C2H2 ), hydrogen ( H2! 
and carbon monoxide ( CO ). The choice of acetylene was based on species 
profiles given for a methane air diffusion flame[12]. Of the intermediate 
hydrocarbons, acetylene was measured to have the highest concentration. 


The general cracking reaction is 


(a oO Pare fy 


th mae A AN Cee i rs Vg es 


with fuel molecules, enthalpies of reaction, and stoichiometric coefficients 
given in table 3. Note that the reaction is endothermic when Qc, is greater 
than zero. The reaction at the flame is then given as 


(bq 
¥ bow Hs ry lt less COne trainee wacianee BO + ae 
C2H, 


This is a very simple approximation to the flame chemistry and will 
probably be improved in the future. But even though it is simple, it will 


be shown to give good results. 


20 


TABLE 2 
Constants for One Step Reaction 


fuel molecule Hy) XK) keo2 420 © Qn) J/k mol) 


cyclohexane (C,H)? l 9 6 6 3.689e+09 

n-Heptane C7H i, l 11 7 8 4.501e+09 

Iso-Octane CsHis 1 12d iS 9 5.100e+09 

Ethanol C2HgO l 3 2 3 1.235e+09 
TABLE 3 


Constants for the Decomposition of Fuel into Acetylene 


fuel molecule Ww Jeou2 4 Ko 3 QOcp VJ/kmol) 


cyclohexane (Hj? l 3. 3 ON 8.03e+08 
n-Heptane CrH 1, | 3.5. <u RO 9.82e+08 
Iso-Octane CgHis l 4 OAS 0 1.13e+09 
Ethanol - C2H,g0 l 0.5 20 l 2.98e+08 


The temperature at which cracking occurs determines the location 
of the cracking sheet. Based on species and temperature data for a 
Wolfhard-Parker slot burner [12], the cracking temperature was chosen to 
be 1350 K. This temperature corresponds to the location of the peak 


concentration of acetylene. 


The ambient parameters are given in table 4. All values are taken 
for air at the ambient temperature except the specific heat and the 
Prandtl number. The specific heat was chosen so that the theoretical 
temperatures do not exceed the adiabatic flame temperature for the fuel. 
_ As areference, this value of specific heat corresponds to the specific heat 


of air at about 2200 K. The presence of species related to combustion 


Z| 


might explain the high value of the specific heat. The Prandtl number was 
increased to 0.73 to account for products in system. The dimensionless 
numbers which control the diffusion flame are given in table 5. The 
definitions for the dimensionless numbers are given in table 1. For the 
mass transfer number, the latent heat includes the heat required to bring 
the liquid fuel from the initial reservoir temperature to the boiling 
temperature. The data for the nondimensional numbers were taken from 


the CRC Handbook of Chemistry and Physics [13]. 


The flame model was run with and without cracking. For the case 


without cracking, Q was set equal to zero and Nc, was set equal to the 


flame enthalpy. The resulting solution 1s equivalent to the case without 


cracking. The temperature dependent density viscosity was used for all 


runs. E 
Table 4 
Values of Ambient Parameters 
Pr 0.73 
Xs 0.232 
it 288 K 
V 146x110 m/sé 
Cp 1405 J/( kg Kk) 
Table 5 | 
Dimensionless Flame Controlling Numbers 
fuel eT B pee fe had gy wast Gry 
cyclohexane 14 7rackin24d 020.32 onc 0-50 0.22 1.08 x10!9 
n-heptane [5H 4.50 0.040 0.50 0.22 1.37110!9 
isO-octane pul 4.95 0.041 0.50 0.22 1.37 x10! 


ethanol 8.98 3.04 0.030 0.50 0.24 1.03 x10! 


ay 


It is useful to collapse the experimental measurements taken at 
different locations downstream onto a single graph through the use of a 
similarity variable. The analysis predicts that the boundary layer will 
grow proportional to the fourth root of downstream location. The 
similarity variable used in the analysis includes the Horwarth transform 
which distorts the shape of the profiles and is unsuitable for use with 
experimental! data. A similarity variable which can be determined from 


experimental measurements alone is defined‘as 


K Eyl int (70) 


Ox» FF x Se e 


and is consistent with the analysis. The analysis predicts that the 
streamwise velocity will grow proportional to the square root of 
downstream location. The similar velocity is defined as 


ad, 
Vx : (7/7) 


= 
— 


To represent the pyrolysis rate at the wall a nondimensional mass flux is 
introduced. From the similar stream function and the nondimensional 
normal velocity, the nondimensiona!l mass flux is defined as 

° 7 ree Le pte bn (é) 


aa re: apa e ae = — CAL 
) aan flo) ) 


hy 
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B. Experimental Results 


The results for temperature and streamwise velocity fields are 
given in figures 2 through 9 using the flame controlling numbers given in 
tables 5 and the ambient properties in table 4. The experimental data 
were taken from reference [2]. The curves labeled one step represent the 
solution to the problem without thermal! cracking. The curves labeled two 
step represent the solution of the problem including the effects of 
cracking. The data points labeled in centimeters indicate at what location 


downstream the data was measured at. 


In general, the inclusion of thermal cracking significantly improves 
the agreement between theoretical and experimental results. The , 
temperature profiles with the effect of cracking exhibit remarkable 
agreement with experimental. For ethanol, n-heptane, and 1so-octane, the 
predicted flame location and the experimental flame location agree to 
Within 2%. For cyclohexane the predicted flame location was off by about 
10%. For the solutions without cracking, the flame locations were off bv at 
least 25%. The shape of the theoretical temperature curves with thermal 
cracking gives a much better fit to the experimental results then the 
temperature curves without cracking. The improvement is particularly 
evident in the temperatures next to the wall. The analytical mode! 
predicts that the cracking sheet will appear as a kink in the temperature 
profile. The experimental temperatures do not show any sort of kink. 
This sugests that the cracking reaction must occur over a wide range of 
temperatures and should not, in general, be confined to a sheet. The 


advantage of the cracking sheet is that it makes the solution far more 
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tractable and is therefore, a useful idealization. The flame is well 
represented by a sheet. The experimental temperatures show a distinct 


peak at the flame. 


The density of the gases between the wall and the flame have a - 
strong inf luence on the flame location due to the appearance of density in 
the Horwarth transformation. The lower temperatures between the wal! 
and the flame cause the flame to move toward the wall! due to the 
increase in density of the gases in this region. This is the main reason 
why the location of the flame moves closer to the wall with the addition of 


cracking. 


The lower temperatures between the flame and the wall have a 
strong influence on the velocity profiles since the flow is buoyancy driven. 
The inclusion of cracking reduces the peak velocity by about 15 %. 
Agreement is improved with velocities on the edge of the boundary layer 
but theoretical velocities closer to the wall are consistently lower than the 
experimental velocities. The experimental velocities do not collapse as 
well as the experimental temperatures. The velocities at two centimeters 
downstream are much higher than expected. It would be useful to have 
more velocity measurements further downstream to establish some sort 


of trend. 


The variable properties allow the mode! to accurately predict the 
flame location. If the density viscosity product (pi) were taken as a 
constant and the ambient values were used, the flame location would be. 
moved outward by about 25%. This would destroy the agreement 


between the theoretical profiles and the experimental profiles. 


25 
C. Pyrolysis Rates 


The agreement between experimental and theoretical temperatures 
near the wall should improve the prediction of the pyrolysis rate at the 
wall since the wall mass flux is determined by the heat transfer to the 
wall. The dimensionless mass flux at the wall 1s graphed for each of the 
fuels in figure 10 comparing the mass flux with and without the cracking 
reaction. The endother mic cracking reaction reduces the pyrolysis rate at 
the wall by about 20%. This result is observed in the temperature plots as 


a decrease in the slope of the temperature profile at the wall. 


cyclohexane n-heptane  iso-octane ethanol 
Fig. 10 


D. Parametric Studies 


Two parametric studies are done using the model to explore how 
the cracking parameters affect the flame. The dimensionless cracking 


enthalpy determines where the cracking reaction will occur. Figure 1] 


shows the results for four different her with all of the other flame 


controlling nu mbers held constant. Decreasing h,, moves the flame closer 


to the wall due to the decrease in temperatures between the flame and 


the wall. In going from Ner= 0.83 to Hee - 0.31 the wall pyrolysis rate 
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drops more than 30 %. This is indicated by the decrease in the slope of 


the temperature profile at the wall. 


Figure 12 shows the effect of changing Q on the temperature 
profiles. Increasing Q moves the flame closer to the wall and causes the 
kink in the temperature profile to become more pronounced. As Q 
increases, the heat flux to the cracking sheet must increase. This causes 
the temperature profile to become steeper on the flame side of the 


cracking sheet. 


_ Conclusion 


Endother mic cracking of fuel and temperature dependent density- 
viscosity product were added to the free convective boundary layer 
diffusion flame solution presented by Pagnil2]. With the addition of 
thermal cracking, the agreement between experimental and theoretical! 
and experimental results are greatly improved. The pyrolysis rate at the 
wall is reduced due the existence of the endothermic peacrine sheet. The 
use of temperature dependent density-viscosity product enables the 
flame model to accurately predict the location of flame. Future work 
would include adding radiation to the mode! or applying the mode! to 


other types of diffusion flames. 


Nomenclature 
ital Letters 

B Spaldings mass transfer number, eqn. (36) 
Bc; Cracking mass transfer number, eqn. (39) 
Berf Effective B number which accounts for cracking, eqn. (40) 
D Diffusivity 
Fy Uncracked fuel 
F2 Cracked fuel molecule 
Gr; Grashof number evaluated atx -1m 
J Solution to energy-species equation 
ee Value of J at crack sheet, eqn. (43) 
Jey Value of J at flame sheet, eqn. (45) 
u Latent heat of vaporization for pyrolyzed material 
Le Lewis number 
Mr Molecular weight of uncracked fuel 
0 Oxygen 0» 
P General term for products 
Pr Prandtl number 
0) Heat of reaction ratio, eqn.(34) 
Qn Heat of reaction for one step reaction 
Qer Heat of reaction for cracking reaction 
Qerfi Heat of reaction for combustion of cracked fuel, eqn. (4) 


T Temperature 


| 
a 
4 


28 
Temperature of fuel reservoir 
Wall temperature 
Temperature at which cracking occurs 


Ambient temperature 
Nondimensional streamwise velocity 
Nondimensional transverse velocity 
Mass fraction of uncracked fuel 


Mass fraction of cracked fuel 


Mass fraction of cracked fuel at wall . 


Mass fraction of uncracked fuel at wall 
Mass fraction of fuel in pyrolyzed material from wall 
Mass fraction of oxygen 


Ambient mass fraction of oxygen 


Small Letters 
Specific heat 
Similar part of stream function 
Enthalpy 
Dimensionless enthalpy, eqn.(34) 
Dimensionless enthalpy at crack sheet, eqn.(44) 
Dimensionless enthalpy at wall, eqn.(36) 
Arbitrary length, set equal tol m 


Mass flux at wall 


Mass generation term 

Dimensionless mass flux at wall, eqn. (73) 
Nor mal to wall 

Pressure 

Streamwise velocity 

Transverse velocity 

Evaluated at wall 

Streamwise coordinate 

Transverse coordinate 

Dimensionless uncracked fuel, eqn.(23A) 
Dimensionless cracked fuel, eqn.(23B) 
Dimensionless oxygen, eqn.(23C) 


Dimensionless cracked fuel at cracking sheet, eqn.(44) 


bols 
Temperature dependent density viscosity function, eqn.(49) 
Dimensionless fuel in transferred gas, eqn.(37) 
Dimensionless uncracked fuel equivalence ratio, eqn.(24) 
Dimensionless cracked fuel equivalence ratio, eqn.(25) 
Stoichiometric coefficient 
Density 
Dynamic viscosity 
Kinematic viscosity at ambient 


Ther mal condictivity 
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anny 


Ce a 


CR 


nit 
f2 


Nondimensional streamwise coordinate, eqn.(52A) 
Nondimensional transverse coordinate, eqn.\52B) 
Nondimensional stream function, eqn.(57) 


Nondimensional similarity variable, eqn.(58) 


Subscripts 
Evaluated at cracking sheet 
Evaluated at flame 


Uncracked fuel 


‘Uncracked fuel 


Cracked fuel 
Green | 
products 
Evaluated at wall 
Evaluated at infinity 
Superscripts 


Per unit time 


Per unit length or derivative with respect to the independent 


variable 
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Figure }. 


Figure 2 


Figure 3 


Figure 4 
Figure 5 
Figure 6 


Figure 7 


Schematic diagram of free convective boundary layer diffusion 
flame showing internal parts and coordinate system. 


Temperature results for cyclohexane. Experimental 
temperature data is compared with theoretical results for the 
model with and without cracking. The curve labeled ‘one 
step” indicates the solution without cracking. The curve 
labeled ‘two step’ indicates the solution including the effects 
of cracking. Experimental! data is taken from three locations 
downstream. The use of a similarity variable, eqn.(73), allow 
the data to be plotted on a single graph. The data labeled 2 cm 
indicates that the experimental data was measured 2cm 
downstream from the leading edge of the flame. The points 
labeled 4 cm and.6 cm are experimental] date taken at four 
and six centimeters downstream. 


Streamwise velocity results for cyclohexane. Experimental 
streamwise velocity data is compared with theoretical results 
for the mode! with and without cracking. The curve labeled 
“one step” indicates the solution without cracking. The curve 
labeled ‘two step” indicates the solution includes the effects of 
cracking. Experimental data is taken from two locations 
downstream. The use of a similarity variable, eqn.(73), and a 
similar velocity, eqn.(74), allow the data to be plotted on a 
single graph. The points labeled 2 cm and 4 cm are 
experimental date taken at two and four centimeters 
downstream from leading edge of flame. 


Temperature results for n-heptane. Symbols and curves are 
defined as in figure 2. 


Streamwise velocity results for n-heptane. Symbols and 
curves are defined as in figure 3. 


Temperature results for iso-octane. Symbols and curves are 
defined as in figure 2. 


Streamwise velocity results for iso-octane. Symbols and 
curves are defined as in figure 3. 


Figure 8 
Figure 9 


Figure 10 


Figure 11 


Figure 12 
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Temperature results for ethanol. Symbols and curves are 
defined as in figure 2. 


Streamwise velocity results for ethanol. Symbols and curves 
are defined as in figure 3. 


The value of the dimensionless mass flux at the wall is 
graphed for each of the fuels for the flame with and without 
thermal cracking. 


Influence of h,, on temperature profiles. The cracking 
temperatures are as follows: h,, 70.83 , Ter = 2050 ‘Hep songs 
, Ter = 1700; he, =0.50 , Ter = 1350: h,, 70.31 . Ter = 950. 
The other flame controlling numbers were held constant. f FT" 
14.7, B= 5.24, Ay = 0.03, J = 0.22, Gry = 1.08 x 1010. 


Influence of Qon temperature profiles.The other flame 
controlling numbers were held constant. Fey = 14.7, B=hiea 
Ay = 0.03, Hcp = 0.50, Gr; = 1.08'x 1019. 
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